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Abstract 

We construct steady states of the Euler-Poisson system with a 
barotropic equation of state as minimizers of a suitably defined en- 
ergy functional. Their minimizing property implies the non-linear 
stability of such states against general, i.e., not necessarily spheri- 
cally symmetric perturbations. The mathematical approach is based 
on previous stability results for the Vlasov-Poisson system by Y. Guo 
and the author, exploiting the energy- Casimir technique. The analysis 
is conditional in the sense that it assumes the existence of solutions 
to the initial value problem for the Euler-Poisson system which pre- 
serve mass and energy. The relation between the Euler-Poisson and 
the Vlasov-Poisson system in this context is also explored. 

1 Introduction 

Consider a self-gravitating fluid ball in IR^ where the fluid has mass density 
p = p{t,x) >0 and velocity field u = u{t,x) & JR^; t>0 denotes time and xE 
JR^ position. In the simplest case such a system obeys the Euler-Poisson 
equations 

dtp + V-{pu) = 0, (1.1) 
pdtU + p{u-\/)u = -\/p-pVV, (1.2) 
A\/ = 47rp, lim V{t,x) = 0, (1.3) 
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where V = V(t,x) denotes the self-consistent gravitational potential and p = 
p{t,x) is the pressure in the fluid. To close this system an equation of state 
needs to be supplied which relates the pressure and the density: 

p = P{p). (1.4) 

An equation of state where the pressure does not depend on the temper- 
ature or specific entropy is often referred to as barotropic. More sophisti- 
cated models for a star would contain equations for the temperature and for 
the thermonuclear processes within the star, but in astrophysics the present 
model is often used as a simple description of a gaseous star, in particular 
if — as in the present note — the issue is the stabihty of steady states since this 
problem is non-trivial already for this simple model. 

To approach this problem we define the energy of a state {p,u) as 

n{p,u):=l f \u\^pdx+ f ^{p)dx-\ f ^^f^^dydx. (1.5) 

Here $ is defined in such a way that for p > 0, 

P\p) = pV{p). (1.6) 

Under appropriate assumptions on $ or P respectively, which include — but 
are much more general than — equations of state of the form P{p) = cp"' with 
7 >4/3, this functional has a minimizer under the mass constraint 

j pdx = M 

with M > prescribed. This follows from a result of the author established in 
PB| . Such a minimizer is a spherically symmetric steady state of the Euler- 
Poisson system with vanishing velocity field and as such is referred to as a 
static solution of the system, as opposed to the more general concept of a sta- 
tionary solution, which is also time-independent but may have non-vanishing 
velocity field. From the fact that the steady state minimizes the energy one 
can deduce a rigorous, non-linear stability result. Since in order to solve the 
variational problem stated above we need to make no symmetry assumption 
we obtain stability against general, i.e., not necessarily spherically symmetric 
perturbations, even though the steady state obtained as an energy minimizer 
must a posteriori be spherically symmetric. An essential point for the stabil- 
ity analysis is that both the energy defined in ( [1.5|) and the total mass J pdx 
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are conserved along solutions of the time dependent problem. Since the ini- 
tial value problem for the Euler-Poisson system is far from being completely 
understood and one has to expect the occurrence of shocks or gravitational 
collapse our stability result is conditional in the following sense: As long as 
a solution to the initial value problem exists and preserves mass and energy 
it satisfies the stability estimate. 

The paper proceeds as follows: In the next section we study the varia- 
tional problem stated above. Then we show that a minimizer is a spherically 
symmetric, static solution. In Section 3 we carry out the stability analy- 
sis. In a final section we discuss the relation of the present result to those 
obtained by Y. Guo and the author for the Vlasov-Poisson system 

dJ + vVJ-VV-VJ- = 0, 

AV = 4'Kp, lim V{t,x) = 0, 

>oo 

p{t,x) = / f{t,x,v)dv. 



Here the dynamic variable is the number density in phase space, / = f{t,x,v), 
of an ensemble of massive particles with spatial density p, interacting by the 
gravitational potential V which the particles create collectively; v e IR^ de- 
notes the momentum or velocity coordinate in phase space. In this model 
collisions among the particles are neglected which is appropriate if one is de- 
scribing a galaxy or a globular cluster. In |^ ||, |^, |10|, ^ a variational 



approach to the question of the stability of steady states of this system was 
developed, which relied on the minimization of appropriately defined energy- 
Casimir functionals. In p3| the author observed that these functionals, which 



act on phase space densities f = f{x,v), can in a natural way be reduced to 
functionals acting on spatial densities p = p{x), and he studied the variational 
problem for these reduced functionals, without realizing that this yields the 
stability result for the Euler-Poisson problem which we discuss here. Inci- 
dentally, the close relation between the stability of barotropic stars and of 
stellar systems was already observed on a formal level in the astrophysics 
literature and is often referred to as Antonov's First Law, cf. 5.2 (b)]. 

We conclude this introduction with some additional references to the lit- 
erature. A lot of background on the physics of stars can be found in p and 
|13[. An excellent and broad overview of mathematical results for hydro- 



dynamical models of gaseous stars is given in M which also contains many 
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further references. Our variational approach is related to the concentration- 
compactness principle due to P.-L. Lions [|17|. The use of energy- Casimir 
functionals for questions of stability was discussed in a very broad context in 
1l2| . As aimed more specifically to the existence and stability properties of 



steady states for gaseous stars we mention M, [TBI and the references there. 



2 Energy minimizers and steady states 

If we minimize the energy Ti as defined in ( [L.5|) over all "relevant" states 
{p,u) then a minimizer will clearly satisfy u = 0. Thus we study the following 
variational problem: For a density p we define the induced potential 

J \x — y\ 

and the potential energy 

^pot(p) ■■=-111 ^^0^'^y^^=l J pVpdx = -^ J \VVp\'dx. (2.1) 
We want to minimize the functional 

Hrip) := I ^pix))dx + E^M (2.2) 

over the constraint set 

^m:={pGL1(1R^)|p>0, y'$(p)<oo, Jp = m} (2.3) 

for M > prescribed and $ satisfying the following 
Assumptions on $: $ G C^([0,cx)[), $(0) = = $'(0), and 

($1) $ is strictly convex. 

($2) $(p) > Cpi+V" for p > large, with < ra< 3, 
($3) $(p) < Cpi+i/"' for p > small, with 0<n'<3. 

Note that pe nL^/^{lR^) and hence 1/peLP(lR3), VVpeLi(]R^) for 
PEJ-'m and any pG]3,12] and gG]3/2,12/5] so that the potential energy is 
well defined. The following result was proved earlier by the author, cf. ||2^, 
Thm. 3.1]; it is included here for easier reference. 
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and 



Theorem 1 The functional Tir is hounded from below on with 
inf jF^^ < . Let {pi)GJ^M be a minimizing sequence of Tir- Then there 
exists a sequence of shift vectors (oj) ClR^ and a subsequence, again denoted 
by [pi), such that for any e>0 there exists R>0 with 

[ pi{x)dx>M-e, ieIN, 
Tpi: = pi{- + ai) ^ Po weakly in L^^^^'^{JR^), i-^oo, 
I Po>M-e. 

JBr 

Finally, 

VVtp, strongly in L^(IR^), i—^oo, 

where Vq = Vp^, and PqEJ-'m is a minimizer ofHr- 

Here and in the foUowing we denote for R>0, 

BR:={xeJR^\\x\<R}. 

Remark 1. As noted in the introduction, the function $ appearing in 
the definition of the energy and the function P defining the equation of state 
(|1.4| ) have to satisfy the relation ( |1.6|) . If we think of P as the basic, given 



quantity then we define 

CP rP'ir) 



$(p):= / / ^-Idrda, p>0, 

Jo Jo T 



and this function satisfies the assumptions ($1), ($2), ($3) provided P'{t)/t 
is integrable on ]0,cxo[ and the assumptions 

(PI) P'>0, 

(P2) P'{t)>Ct^/'^, r>0 large, with 0<n< 3, 
(P3) P\t)<Ct^I'^' , r>0 small, with 0<n'<3 
hold. This is the case in particular for polytropic equations of state 

P{p)=cp^ 
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with c>0 and 7>4/3, in which case 



7-1 

Equations of state where the dependence of the pressure on the density is 
different for large and for small values of the density are motivated from a 
physics point of view. 

Remark 2. The need for the spatial shifts in Theorem |l] is not technical 
since without them the assertion of the theorem is false: Given a minimizer 
Po the fact that the functional Tir is translation invariant implies that by 
shifting Po off to infinity one obtains a minimizing sequence which converges 
weakly to zero and not to a minimizer. 

Now that a minimizer is obtained we show that it is a steady state of 
the Euler-Poisson system; in order to avoid technical discussions concerning 
regularity let us now assume in addition to the above that $gC^(]0,oo[) 
which means P G C^(]0,oo[). 

Theorem 2 Suppose that po G is a minimizer of the functional Tir with 
induced potential Vq. Then there exists a Lagrange multiplier Eq<0 such 
that 

({^')-\Eo-Vo),Vo<Eo, 
,Vo>Eo. 



Po- 



The functions po and Vq are spherically symmetric with respect to some point 
in IR^, as functions of the radial variable po is decreasing and Vq increas- 
ing, poGC(IR^) with finite mass M, VoGC^(IR^) with lim|^|_>ooVo(a;) = 0, 
and (po,Mo = 0) is a time-independent solution of the Euler-Poisson system 
with equation of state 

P{p)= ra^"{a)da. 
Jo 

If^'{p) <Cp^/'^' for p>0 small then Eq<0 and po has compact support. 

The assumptions on $ imply that $' is invertible on [0,oo[. The growth 
condition on $' implies the condition ($3) but is not equivalent to it. 
Proof of Theorem |2|. We start by deriving the Euler-Lagrange equation for 
the variational problem. Let po G J-'m be a minimizer with induced potential 
Vq. For e > define 

^,:={a;GlR^|e<po(x)<l/e}; 
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think of po a pointwise defined representative of the minimizer. For a 
test function w G L°°(IR^) which has compact support and is non-negative on 
\ Se define for r > small, 

Jwdy 

pT- = Po + rw-T——ls,, 

where 1^^ denotes the indicator function of the set S^- Then Pr>0 and 
J pr = M so that pr G for t > small. Since po is a minimizer of Tir, 

< nr{Pr)-nr{po) =rj ($'(po) ^^^^ " ^^^s}j dx + o{T). 

Hence the coefficient of r in this estimate must be non-negative, which we 
can rewrite in the form 



fWipo) + Vo \-([ i^'ipo) + Vo)dy 

J L vol be \-Js^ 



wdx>0. 



This holds for all test functions w as specified above, and hence $'(po) + Vq = 
E, on S, and $'(po) + Vo>E, on IR^ \ for all e > small enough. Here 
is some constant which by the first relation must be independent of e, and 
taking e — proves the relation between po and Vq. 

The symmetry assertion follows by a rearrangement argument: Let pg 
denote the symmetric decreasing rearrangement of po. Then 



<l>(po)=/$(Po), ^pot(Po)>i^pot(pS), 



and Pq gTm, cf. |T^, Thms. 3.7, 3.9]. Since po is a minimizer of Hr equality 
must hold in the estimate for the potential energy which implies that po must 
be spherically symmetric with respect to some point in IR^, and decreasing 
as a function of the radial variable. Alternatively, one could use Thm. 4] 
to conclude the spherical symmetry. 

We now continue essentially as in the proof of [p2| , Thm. 3]. First we 
observe that by assumptions ($1) and ($2) and the mean value theorem 

$-(p)>$V) = ^(^^-^(°^>CpV" 

p-0 

for all p large, with some intermediate value < a < p. Hence 

Po{x)<C{l + iEo-Voix)r) 
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and since Vq e L'^'^ (IR^) and n<3 we have poEL^nL'^iM^). For any R>1 
we split the convohition integral for Vq according to \x — y\ <1/ R, 1/R< 
\x — y\< R, and \x — y\>R to obtain 

/ /-V-R , \^^'^ r M 
-Vo{x)<C\\po\\A r^-^/^dr +R p^{^y)dy + —,\x\>R. 

\Jq J J\y\>\x\-R R 

This implies that VoeL°°(]R^) with Vo(a;)— >-0, |a;|— >-oo. The asserted regu- 
larity of po and Vq now follows from the relation between these two quantities 
and Sobolev's embedding theorem. 

The limiting behavior of Vq and the relation between po and Vq implies 
that Eq < 0, since otherwise Pq{x) > {^')~^{Eq/2) > for large which would 
contradict the integrability of po- 

Next we check that (po,Mo = 0) solves the stationary Euler-Poisson system 
with the asserted equation of state. Since by construction Vq is the potential 
induced by po all that needs to be checked is the Euler equation which takes 
the form 

Vpq + PqVVq = 0. 

Now ^'{pq) = Eq — Vq on supppo. Taking the gradient of this relation proves 
the assertion, provided P' = $"p. 

It remains to show that under the additional growth condition on the 
density is compactly supported and the cut-off energy Eq strictly negative. 
Since Vq increases to zero at spatial infinity the relation between po and Vq 
implies that po has compact support provided Eq < 0. If we assume that 
£"0 = then the additional assumption on $' implies that 

PQ{x)>C{-VQ{x)f 

for \x\ large, i.e., po small. But a simple expansion of the Green's function in 
the convolution formula for Vq and the fact that Jpo — M show that —Vo{x) > 
M/(3\x\) for l^l large. Inserting this into the estimate for po from below 
yields a contradiction to the integrability of po. Hence £^0 < 0, and the proof 
is complete. □ 

In view of the stability result discussed in the next section it would be 
desirable to know that for fixed M there is up to spatial shifts a unique min- 
imizer of Hr in J^m or at least that the minimizers are isolated up to spatial 
shifts. Numerical evidence obtained by solving the spherically symmetric 
Poisson problem AVo = 47rpo = 47r(<l>')~^(£o~H) seems to indicate that the 
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minimizers are in general not unique but are isolated. An example where the 
minimizer is unique is provided by the polytropic equation of state: 

Remark 3. Let P{p) = cp^^ i.e., ^{p) = ^p^ with 7>4/3. Then for 
every M > there exists up to spatial shifts exactly one minimizer po ^ 
of Tir- This follows from the fact that Eq — Vq solves the Emden- Fowler 
equation 

l-^r''z')' = -czl, r>0, 

where 1 + 1/?t, = 7 and z+ denotes the positive part of z. All solutions of this 
equation which are regular at the center are related by a scaling transfor- 
mation, and prescribing the mass fixes the corresponding scaling parameter. 
For more details cf. |2^, Thm. 3 (b)]. 



3 The stability analysis 

We start by expanding the total energy Ti about a minimizer po- For P ^ 
and u such that / |Mpp < oo we have 

n{p,u)-n{po.Q) = d{p,po)-^\\vVp-vv4l + ]^j H'^pdx (s.i) 



where 



d{p,po):= / $(p)-$(po) + (V^o-^o)(p-Po) 



dx: 



we are allowed to insert the term Eo^p — po) into d since J p = J Po- The fact 
that $ is convex and the relation between po and Vq established in Theorem ^ 
imply that for all pG J-m, 

d{p,po)>0. 

Moreover, d{p,po) = only if p = po, and if one assumes in addition that $" 
is bounded away from zero then 

d{p,po)>C\\p-po\\l, peJ^M- 



What is irritating about is the minus sign in front of the difference term 
for the gravitational fields, but this term converges to zero along minimizing 
sequences, cf. Theorem ||. 

Our stability analysis relies on the fact that the mass and the total energy 
Ti. defined in ( |1.5| ) are conserved along solutions of the time-dependent Euler- 
Poisson system. Formally, conservation of mass follows from the continuity 
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equation ( p. . 1| ) , and conservation of energy follows by a straight forward com- 
putation from all three equations ( p. .21) , ( pTS]) together, provided the 
relation ( |1.6| ) between the equation of state p = P{p) and the function $ 
which appears in the energy holds. However, from a rigorous mathematical 
point of view insisting on these conservation laws amounts to an assump- 
tion on the solutions of the initial value problem. This is reflected in the 
statement of our stability result. 

Theorem 3 Let pqEJ-'m be a minimizer of Ti,- with induced potential Vq, 
and assume that it is unique up to spatial shifts. Then for every e > there 
is a 5>0 such that the following holds: For every solution ty-^ {p(t),u{t)) of 
the Euler-Poisson system with p(0) G Tm which preserves mass and energy 
the initial estimate 

^(p(0),Po) + ^ II Vl^p(o) - VK)||2 + ^ / 1^(0) rp(0) dx < S 

implies that for every t>0 for which the solution still exists there is a shift 
vector a G IR^ such that 

d(p(t),TVo) + ^ II VV;(t) - Wtvo Wl + lJ Ht)\'p{t) dx < e; 
here T"-p\ = p{- +a). 

Proof. Assume the assertion were false. Then there exist e > 0, t„ > 0, and 
initial data p„(0) eJ-'m and u„(0) such that for all nG IN, 

diPniO),po) + ^11 V\/p„(o) - V^olls + ^/ |wn(0) I V(0) dx<^ (3.2) 
but for any a G IR^, 

rf(p„(t„),T>0) + ^ II VV;„(i„) - Wrap, \\1 + \J \Un{tn)\'Pn{tn) dx > t. (3.3) 

By (§3) and 

Jim 7i:(p„(o),u„(o)) =7i:(po,o) =7i:.(po). 

Since by assumption the solution 1 1— > {pn(t) ,u„it)) conserves energy, 

lim sup K(Pn (in)) < lim 'H{pn{tn),Un{tn)) ='Hripo), 
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and since it conserves mass, (p„,(t„)) dJ^M is a minimizing sequence for Ti^- 
By Theorem |l] there exists a sequence (a„) C IR^ such that up to a subse- 
quence, 

l|V\^p„(t„)-V\/T^„pJ|2--0; (3.4) 

at this point we used the uniqueness of the minimizer. Note also that for any 
p E Tm and a G IR^, 



|VVV.,-VK 



P0II2- 



Since hm„^oo^(Pn(i^n),Mn(i^n)) =^r(Po) =^r(7'""Po) we couclude by 
and ( |3.1| ) that 

d{pn{tn),T"-"pQ) + ]^j \Un{tn)\^Pn{tn)dx-*Q, U^OO, 



a contradiction to (|3.3|). □ 

Remark 4. Without the uniqueness assumption for the minimizer we 
obtain a stabihty resuh of the following type: Let J^m C. Tm denote the set 
of all minimizers of in Tm- Then for every e > there is a 5 > such that 
for any solution with p(0) G Tm the initial estimate 



inf 



rf(p(0) ,po) + ^ II VV^p(o) - Wo II 2 + ^ / 1 w(0) rp(0) dx 



<6 



implies that 



inf 

Pq<^Mm 



d{pit),po) + ^ II VI^pW -VVj\l + y \uit)\'pit) dx 



<e 



(3.5) 



as long as the solution {p{t),u(t)) exists and preserves mass and energy. The 
proof of this assertion follows exactly the same line of reasoning as the proof 
of Theorem ^ 

Assume now that the minimizer po is not necessarily unique, but isolated 
up to spatial shifts, that is to say, 

60 : = inf { II V - V^p,, II 2 I Po G \ {r>o | a G IR^} } > 0. 

In this case the assertion of the theorem again holds, provided the solution 
of the Euler-Poisson system is sufficiently continuous in t so that it cannot 
jump from one minimizer to the next, more precisely: Let e>0 arbitrary. In 
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order to find the corresponding S we can without loss of generality assume 
that e<(5o/4. Now choose 6>0 according to the first part of the present 
remark so that ( p.5[ ) holds, without loss of generality 6 <e, and let p(0) G J-'m 
and u{0) be such that 

c?(p(0),po) + ^||VV,(o)-VK)||^ + i / \u{0)\'p{0)dx<6. 

OTT Z J 

The required continuity assumption on the corresponding solution is that 

h{t, a) := d(p(t),T>o) + ^ II Vl^p(t) - VVt^p, II2 + ^ / 1«(^) I dx 

is continuous, and indeed, that inf^^-^3h{t,a) is continuous. Now assume 
that there exists t>0 such that 

inf h(t,a) > e. 

Since at time zero the left hand side is less then e there exists some t* >0 
where 

inf h{t*,a) = e. (3.6) 

On the other hand, the first part of the present remark provides some Pq G 
Mm such that 

«n,Po) + ^l|VV;(t.)-Vy,s||^ + ^/ \uif)\'pif)dx<e<^-^. (3.7) 
By ( |3.6|) and ( p.7|) together with the non-negativity of d, 

^l|Vt'„-VV',.,||^<|, 

and by the definition of 5q there must exist some a* G IR^ such that Pq = T""* pq. 
But this means that ( |3.6|) contradicts ( |3.7|) , and the argument is complete. 

Remark 5. The restriction p{0)eJ-'m for the perturbed initial data is 
acceptable from a physics point of view: A small perturbation of a given 
star, say by the gravitational pull of some outside object, results in a per- 
turbed state with the same mass. However, such a perturbation will hardly 
be spherically symmetric so that it is important that no such restriction is 
necessary in our stability result. On the other hand, if we restrict ourselves to 
spherically symmetric perturbations then the spatial shifts in the statement 
of the stability result are no longer necessary. Their need arises from the fact 
that a given steady state can be given a uniform velocity in one direction so 
that it drifts off and the distance of the perturbed state from the original one 
grows linearly in t, no matter how small the initial perturbation was. 
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4 Fluid models versus kinetic models 



Obviously, the hard part in the analysis, if any, is the proof of Theorem]^. The 
essential arguments for that theorem as well as its exploitation for stability 
questions arose from the investigation of the stability of galaxies by Y. Guo 
and the author. In typical galaxies even stars which are spatially close to 
each other can have very different velocities, and hence galaxies are modeled 
not by fluid equations but by kinetic equations, i.e., by the Vlasov-Poisson 
system which was stated in the introduction. On the level of existence and 
regularity of solutions to the initial value problem the Vlasov-Poisson system 
is very different from the Euler-Poisson system: Continuously differentiable 
and compactly supported initial data for / launch continuously differentiable 
solutions which are global in time, cf. [IIU], ITB|, p3, and neither shocks nor a 



gravitational collapse can occur. On the other hand, as to the existence and 
stability of stationary solutions both systems seem to be intimately related, 
and in this last section we want to comment on this relation. 

In [|, 1^, H, I], ^ steady states of the Vlasov-Poisson system were 



obtained as minimizers of an energy- Casimir functional 
nc{f) = JjQ{f{x,v))dvdx + yj\v\'f{x,v)dvdx-yj^4^^^dxdy 



on the constraint set 



Gm ■■= {/ G L\]R') I / > 0, C(/) +Eki„(/) <oo, JJ fdvdx = m}, 

where C(/) and -Ekml/) denote the first and second term in the energy- 
Casimir functional respectively. Such minimizers exist, satisfy assertions 
analogous to Theorems |1| and 0, and are non-linearly stable, provided Q 
satisfies assumptions which are exactly parallel to the ones for $, except that 
the parameters < n,n' < 3 have to be replaced by parameters 0<k,k' < 3/2. 



In it was observed that there is a one-to-one correspondence between 
the minimizers of the energy-Casimir functional Tic on the "kinetic" set Qm 
and the minimizers of the "reduced" functional TCr on the set J-'m- If Po is a 
minimizer of the latter functional with induced potential Vq then 



/o: 



{Q'r\Eo-E) ,E<Eo, 
,E>Eo 
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is a minimizer of Tic, where the particle energy E is defined as 

E = E{x,v) ■.= ^\v\^ + Vo{x). 

On the other hand, if /o is a minimizer of Tic then the induced spatial density 
Po = J fodv minimizes Tir- Of course, in order for this correspondence to be 
true the functions Q and $ have to be in the proper relation: 

$(r)= mfj (^^\v\^g{v) + Q{g{v))^dv 

where for r > 0, 

a.:={^7G^'(IR')|^?>0, J i^\\v\'g{v)+Q{g{v))^ dv<^, J g{v)dv = ry 
With this relation it follows that for every / G Qm, 

nc{f)>nriPf), 

and if / = /o is a minimizer of Tic over J-'m then equality holds; one first 
minimizes over all / which give the same spatial density p and then minimizes 
over the latter. Indeed, there is a more explicit relation between Q and $ 
via their Legendre transforms; for a function /i: H— >] — cxd,cxo] we denote by 

h*{X) := sup(Ar — h{r)) 
its Legendre transform. If we extend $ and Q to ] — oo,0[ by +00 then 



*-(,./.-(.-1m 



dv 



for A G IR, and it can be seen that the assumptions on Q translate into the 
corresponding ones for $ where the relevant exponents are related by n = 
k + 3/2. These assertions are estabhshed in cf. also [p7| ]. 

If one is interested only in the existence of spherically symmetric steady 
states and not their stability the relation becomes even more direct: Suppose 
/o is such a steady state of the Vlasov-Poisson system and define its radial 
and tangential pressures as functions of the radial variable r= \x\ by 



Po(r) = y"(^^^ fo{x,v)dv, Pair) 



2. 



/ 


XXV 




r 



2 



fo{x,v)dv. 



14 



Then a simple computation using the Vlasov equation imphes the relation 



Po = l{Po-Po)-poV^, (4.1) 

the Tolman-Oppenheimer-Volkov equation. If we make the isotropic ansatz 

foix,v)=^{E{x,v)) (4.2) 

with > prescribed and the particle energy E given in terms of the potential 
Vo as above then the pressure is isotropic, Po=po, and ( [4.1D reduces to the 
static Euler equation p'q + pqVq = 0. Moreover, 

roo 

Po(r) = 2'/\ / <l,iE)iE-Voir)y/'dE=:g4Voir)), 

J Voir) 
97/2 poo 

Po{r) = -^rr / <j){E){E -Vo{r)f/' dE = :h^{Vo{r)). 

3 J Voir) 



Hence the ansatz (|4.2| ) reduces the stationary Vlasov-Poisson system to the 
semilinear Poisson equation 

AVo = A7rg^{Vo). (4.3) 

A solution of the latter automatically induces a solution of the static Euler- 
Poisson system which equation of state PQ = h^og'^^{pQ)] both g^f, and 



are invertible on their support under mild assumptions on 0, cf. On 
the other hand, starting from the static, spherically symmetric Euler-Posson 
system one can integrate the Euler equation using the equation of state and 
express po as a function of the potential Vq so that again the problem is re- 
duced to a semilinear Poisson equation, which is precisely ( [4 .31 ) if the equation 
of state has the form Po = h^og'^^(pQ). 

However, these parallels between static solutions of the two systems break 
down at several points: First of all, static solutions of the Euler-Poisson sys- 
tem, that is, steady states with vanishing velocity field, must be spherically 



symmetric, cf. [1^1, but for the Vlasov-Poisson system there exist axially 
symmetric steady states which are not spherically symmetric and which have 
non- vanishing velocity field Jvfodv^O,ci. pT| . 



Another point where the two systems differ, now concerning the ques- 
tion of stability, is the following: For the Euler-Poisson system the threshold 
n<3 which corresponds to 7 > 4/3 for the equation of state is sharp, since for 
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the polytropic case with 7 = 4/3 the energy of a steady state is zero and an 
arbitrarily small perturbation of such a state can make the energy positive 
and cause part of the system to travel off to infinity. These assertions are 
shown in p, Thm. 1.3 (iv), Thm. 1.4]. For the Vlasov-Poisson system one 
can go beyond this threshold to obtain stability provided 0</c<7/2 which 
corresponds to3/2<'n, = fc + 3/2<5. This was done in [jlO[ by minimizing the 
energy Ei^i^ + Epot under the mass-Casimir constraint / f + C{f) = M. Mov- 
ing the Casimir functional C from the functional which we minimize into the 
constraint allows for an extension of the stability result which includes all 
the polytropes with index up to and including n = k + 3/2 = 5. The reduction 
mechanism which takes us to the Euler-Poisson system does no longer work 
in this context, / ^{p)dx is not conserved by itself so that we cannot move 
this functional from the minimized functional into the constraint, and all 
this fits with the fact that — as opposed to the Vlasov-Poisson system — the 
threshold < 3 is sharp for the Euler-Poisson system as regards stability. 
To conclude we note that if the variational formulation is chosen as in 



10| then one can prove the stability result in the form stated in Theorem 



without assuming uniqueness or even isolatedness of the minimizer, cf. p6 
Whether this is possible also in the framework of the present note is one of 
the many open problems in this area. 
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